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Wave instabilities and breaking. Free-surface 
capillaries (images at left), at a Froude number FR = 
0.43, exhibit large spatial variations of instantaneous 
velocity and concentrations of positive (gray) and 
negative (white) vorticity beneath their troughs and 
crests. An increase to Fr = 0.56 generates a fully- 
evolved breaker, represented by the instantaneous 
velocity in the laboratory frame (top image), and a 
frame moving at one-third the inflow velocity (midcle 
— Instantaneous vorticity concentrations are 
evident in the separated mixing-layer formed by the 
onset of breaking (bottom image). From the labora- 
tory experiments of Lin and Rockwell (Physics of Flu- 
ids 6(9), p.2877 and JFM 302, p.29). 
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FOCUSED ISSUE ON NONLINEAR OCEAN WAVES 


Guest Editors: Louis Goodman 
Michael Shlesinger 
Thomas Swean 


Most of the surface of the earth is covered by water and 
most of the time this surface is in motion in the form of water 
waves. Surface waves influence the performance of all Naval 
Systems that operate on, through, or near the ocean interface 
and are central to a wide range of geophysical issues as well. 
Nonlinear wave behavior is important to nearly all Navy 
concerns regarding surface waves including remote sensing, 
acoustic propagation, extreme wave loads, reliability and sea- 
worthiness, sediment transport and near-shore operations. Im- 
provements in Naval capabilitics to operate in the ocean 
environment depend significantly on advancing the under- 
standing of nonlinear ocean wave generation and evolution. 
This is the subject of this focused issue of the Naval Research 
Reviews. 

The mathematical theory of surface waves has been ac- 
tively developed since the 1800’s. The theory is complex. 
Fundamentally, this is because the instantaneous shape of the 
interface is both a boundary condition and a solution to a 
nonlinear partial differential equation. Furthermore, the time- 
dependent shape of the ocean surface is affected by wind 
conditions, underwater currents, internal waves and bottom 
topography. Classical approaches have gencrally exploited the 
underlying linearity of the nonlinear system, weak nonlineari- 
ties, and specialized nonlinear solutions such as solitary wave- 
forms. It is now recognized that further advances in the 
fundamental issues of ocean wave dynamics necessitates in- 
vestigation of problems that are often highly nonlincar, and 
even turbulent and chaotic. The significance and challenge 
presented by these nonlincar problems resulted in the devel- 
opment of a five-year program on nonlinear ocean waves 
which began in 1992. Major topics of interest have been wave 
generation and growth, wave-wave interaction, wave-seabed 
interaction and wave breaking. New theoretical approaches 
based on extending and applying the ideas of nonlinear dy- 
namics and chaos have been brought to bear on these problems. 
Emerging data analysis techniques (e.g. wavelets, inverse scat- 
tering transforms, higher order spectral transforms) have been 
further developed and applied to field and laboratory data as 
well as to solutions of nonlinear partial differential equations 
governing nonlinear wave behavior. 

Our issue starts off with a beautiful review by Professor 
Mei on the many time and length scales found in ocean wave 
dynamics. He summarizes some of the recent advances in 
understanding the interactions among different scales and 
stresses the practical implications this new understanding por- 
tends for coastal engineering. Dr. Huangs article tells of the 
importance of laboratory wave tank experiments in identifying 
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and studying many of the controlling mechanisms of wave 
dynamics. He reminds us that most theoretical advances have 
been inspired by direct observation of wave phenomena and 
that cleverly designed laboratory experiments constitute the 
bulk of these observations. 

Our next three articles review important analysis tech- 
niques that have been applied to the ocean surface wave 
problem. Professor Spedding reviews the application of two 
dimensional wavelet transforms for analyzing ocean waves. 
The wavelet has the advantage of providing local information, 
whereas the conventional power spectrum localizes in fre- 
quency but mixes all spatial information. Professors Elgar, 
Herbers and Guza review bispectral techniques to study cor- 
relation and energy transfer in triads of waves and discuss their 
application to near-shore wave dynamics to gain insight into 
the nonlinearities of wave interactions. Professor Osborne 
reviews aspects of inverse scattering theory and shows how 
new developments have resulted in remarkable simplifications 
to the use of inverse scattering transforms for the study of 
nonlinear water waves. He uses this technique to identify 
soliton behavior in shallow water waves. This information is 
not readily available or accessible through other approaches. 
Together these techniques provide the modem researcher with 
a new and more powerful array of tools with which to look at 
ocean wave data from a fresh perspective. 

With the exception of studies directed specifically at the 
generation of waves by wind, the classical theories have gen- 
erally ignored the dynamics of the air and neglect the possibil- 
ity of shear in the water. However, it has become recognized 
that there may be a number of important phenomena where air 
and water shears play major roles. Our final article, contributed 
by Professors Caponi and Saffman, describes recent major 
advances in understanding the effect of wind and water shear 
on surface waves. The interplay of long and short waves with 
wind and shear generates a variety of nonlinear behaviors 
whose understanding is important for interpreting radar im- 
ages of the ocean. 

With each issue of the Naval Research Reviews we high- 
light a person who has contributed substantially and consis- 
tently to the subject at hand. In this issue we honor Professor 
Marshall Tulin with a brief biographical note and a sketch of 
his likeness. 

The articles in this issue of the Naval Research Reviews 
represent but a sampling of this work; however, we hope they 
will provide the reader with a good overview of the state-of- 
the-art understanding of the complex behavior of nonlinear 
waves. It is certain the reader will develop an increased appre- 
ciation of ocean wave dynamics as one of the most exciting 
and sophisticated areas of scientific investigation, and one with 
immense practical importance as well. 





Multiple-scale Dynamics of 
Ocean Surface Waves - 


Physics and Engineering 


Applications 


Chiang C. Mei 


Department of Civil and Environmental Engineering 


Massachusetts Institute of Technology 
Cambridge, MA 


Introduction 


One of the most distinguishing features of ocean surface 
waves is the existence of many contrasting length and time 
scales. If one scale is examined at a time, the physics can be 
easily studied, either by theoretical means or in the laboratory. 
The interactions among different scales are, however, basic to 
the variety and complexity of natural phenomena. In the past 
few decades, giant strides have been made on these interac- 
tions, especially when nonlinearity plays a role. 

In this review we shall summarize some of the interesting 
advances, with particular emphasis on the dynamical effects 
of weak nonlinearity, whose influence is, however, not always 
weak over sufficiently long temporal or spatial scale. Selected 
topics relevant to both oceanographic and engineering interests 
will be discussed, with a view to highlighting not only the 
commonality of scientific ideas and techniques but also to 
sitimulating further interactions between ocean scientists and 
engineers. 


We begin with the interaction between short and long 
waves, a topic of key importance to the development of the sea 
spectrum. A partial impetus for studying this topic stems from 
the need for proper interpretation of satellite images from radar 
signals backscattered from the sea surface. These radar waves 
are usually in the wavelength range from centimeters (X-band) 
to tens of centimeters (L-band), and are too short to penetrate 
much below the sea surface. Nevertheless, revelations of de- 
tailed bathymetric features of the sea bottom by these images 
suggest the possible role of long-scale modulation of short 
surface waves by long waves or currents, which are directly 
affected by bathymetry. 

The physics of short and long waves is relatively easy 
to analyse if one or both are of small amplitude. Much 
insight can be gained by studying either weak long waves 
generated by much stronger short waves, or short waves 
riding on stronger long waves (with currents being a limit- 
ing case). These topics are separately discussed in the sec- 
tions to follow. 
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Figure 1. 


Wave spectra recorded at Barbers Point Harbor, Hawaii, 1989, 
by Hawaii District Office, Corps of Engineers. Records in the 
lower (higher) frequency half of the diagram were taken from 
a station inside (outside) the basin. 
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While the science of ocean waves captivates oceanog- 
raphers and engineers alike, the latter are motivated more 
often by phenomena affecting or affected by man-made 
structures. Practical interests in coastal engineering require 
the understanding of long-period resonance of harbors 
when the incident wind waves have much shorter periods. 
Prediction of slow-drift oscillations of tension-leg plat- 
forms must be based on a sound understanding of long-pe- 
riod waves affected by the interaction of short-period wind 
waves and the floating structure. In the last three sections 
we discuss some nonlinear wave problems arising in 
wave/body interactions. 

The dynamics of ocean surface waves is a broad topic. 
Several important aspects not covered here have fortunately 
been surveyed in scattered articles. For example the dynamics 
of the breaking waves in the surf zone has been reviewed by 
Battjes(1988). Many topics of coastal wave dynamics includ- 
ing the formation of sand ripples are surveyed in Mei and 
Liu(1993). Recent developments on wave transformation from 
deep to shallow water can be found in Liu(1994). For deeper 
and more gencral expositions one may consult Phillips (1978) 
for oceanographic applications, Mei(1989) for coastal and 
ocean engineering, and Debnath (1994) for recent mathemati- 
cal advances. 
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Short and Long Waves 


Long waves generated by short waves 


For moored ships, tethered offshore submercibles, or ten- 
sion-leg platforms, it is important to know when to expect 
resonant excitation by incident seas. Since the natural periods 
of many mooring systems often fall in the range of 1 to 10 min., 
sea waves of similar periods, called infragravity waves for 
their much greater length compared to the usual wind waves, 
are therefore of considerable interest to engineers. How are 
these infragravity waves originated? How do they propagate 
over complex bathymetry or coastlines? Do they affect the 
evolution of bathymetry or coastlines? 

One source of infragravity waves is nonlinearity. In par- 
ticular, through second-order interactions any two waves of 
different frequencies from the energetic part of the wind-wave 
spectrum can force waves at the sum frequency @, ,@, and at 
the difference frequency @,_@,. If the spectrum is narrow- 
banded, the difference frequency corresponds to a very long 
period, which would induce long waves in finite depth. In their 
classic theory on radiation stresses, Longuet-Higgins and Ste- 
wart (1960) have shown that a slowly modulated group of short 
waves must be accompanied by a mean-sea-level setdown, 
whose (negative) amplitude is proportional to the square of the 
local wave amplitude and moves at the group velocity of the 
short waves. The set-down is therefore a long wave bound to 
the envelope of the short waves. 

The radiation stresses are essentially the averaged rate of 
momentum fluxes due to wave fluctuations when the averages 
are taken in depth and over the period of the short waves. The 
basic idea is similar to Reynolds equations by averaging over 
random fluctuations of turbulence. In particular, when the 
short waves are gentle, the equation governing the long-wave 
surface displacement is of the form: 


a rere 
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where S;, are the radiation stress components, being the corre- 


lations of velocity components. In turbulent flows the 
Reynolds stresses are not known a priori, and additional hy- 
potheses are needed to enable the computation of mean-flow 
quantities. Whenever breaking is absent and wave amplitudes 
are small, the radiation stresses can be calculated by using the 
linearized theory of waves. Thus, the long-wave displacement 
must be of second order in wave steepness and can be solved 
as a linear problem with known forcing. 

Among the earliest applications of this theory we may 
include surf beats in one dimension (Longuet-Higgins and 
Stewart, 1960) and the long-shore current in the surf zone, 
where empirical assumptions are added for the breaking waves 
(Bowen, 1969; Longuet-Higgins, 1970a,b). The mean-flow 








Reflection coefficient of waves over sincusoidal bars on a 
horizontal bed. In the experiments by Heathershaw (1982) 
sinusoidal bars of fixed amplitude b and 1m wavelength cover 
the region from x = 0 to 10 m. Incident wavelength is 2 m. Four 
different water depths h were tested. The theory of Mei (1985) 
is shown in solid lines. 








= 7 
b/h = 0.147 
7 


Geter? 





0.127 
a 4204004 


0.10 
eerteeeee”” 





|R (x)| 





TLEST 8 


12 16 








problems are reducible to steady states. Molin (1982), Mei and 
Benmoussa (1984), and Liu (1989) have applied the theory to 
narrow-banded waves refracted by slowly varying submarine 
ridges and canyons. They found that in a shoaling water the 
incident wave group and the long setdown wave, which are 
aligned in constant depth, propagate in different directions 
upon entering the zone of variable depth. Furthermore, there 
is another long wave propagating at the speed of Veh of the 
linearized long waves. The second long wave corresponds to 
the homogeneous solution of ( 2.1) and may be termed the free 
long wave, to be distinguished from the group-bound long 
waves. Free long waves can propagate into regions where there 
are no short waves, whose absence may be the consequence of 
reflection or diffraction by bathymetry or lateral boundaries. 
A number of interesting scenarios may occur when groups of 
short waves propagate from a zone of constant depth to a zone 
of variable depth. For example, if the short waves of a certain 
carrier frequency advance obliquely toward a deep canyon, 
they can be first refracted and then reflected at a caustic so that 
they do not cross the canyon. However, free long waves can 
enter and pass the canyon. In the case of a submarine ridge, the 
short waves pass the ridge but free long waves can be trapped 
on the ridge (Mei and Benmoussa, 1984; Agnon and Mei, 
1988). Short waves incident on a large island should not be 
significant deep inside the shadow, but the long waves can 
creep in there (Zhou and Liu, 1987). If groups of short waves 
enter a harbor, free long waves can be created at the harbor 
entrance and be resonated inside. This phenomenon has been 


experienced in many harbors where tsunamis are of no con- 
cer. A record of this phenomenon is shown in Figure 1, where 
the wave spectra from gages outside and inside the Barbers’ 
Point Harbor, Hawaii, are superimposed. While the records 
taken at a station outside give only the high-frequency peaks, 
the low-frequency peaks are found at a station inside the 
harbor. Theories for such phenomenon are still at their infancy 
(Mei and Agnon, 1989; Wu and Liu, 1990). 

Along some beaches such as the Chesapeake Bay, there 
are many alongshore bars running parallel to the coast. Bragg 
resonance can cause strong reflections of the incident waves if 
the spacing between bars is roughly one-half that of the incom- 
ing sea waves. This was measured in the laboratory by Heath- 
ershaw (1982) and predicted theoretically by Mei (1985). 
Figure 2 shows the spatial variation of the reflection coefficient 
over ten sinusoidal bars. This strong reflection has at least two 
nonlinear effects. One is to induce a circulation in the bottom 
boundary layer, and one is to create long waves, both bound to 
the short wave envelope and free with the speed Veh. The first 
effect is the result of time-averaged Reynolds stresses due to 
friction in the bottom boundary layer. Very close to the seabed 
the induced drift current converges toward the points directly 
beneath the maxima of the wave envelope; near the top of the 
boundary layer, the drift current converges toward the points 
beneath the minima of the wave envelope. Thus, heavy sand 
particles rolling on the seabed are transported toward the 
envelope minima to form or to heighten the crests of sandbars. 
Since these extrema are separated at half of the incident 
wavelength, the wavelength of sand-bars is roughly one half 
of that of the waves. Light particles in suspension are trans- 
ported toward the minima of the wave envelope. With more 
deposition at the crests, parallel bars grow in height and 
increase the reflection by Bragg resonance, which in tum 
further enhances the growth of existing bars and the birth of 
new bars (see Mei, 1989, p. 433). Thus the increase of wave 
reflection and the growth of sandbars are cause and effect of 
each other. As another consequence, sand particles of different 
diameters are segregated so that coarse sands are more abun- 
dant near the bar crests and fine sands near the troughs; this 
phenomenon is called sediment sorting. 

The second nonlinear effect is the result of radiation 
stresses in the inviscid core, and is present when the incident 
waves are narrow-banded. Then scattering of short waves is 
accompanied by the scattering of bound long waves and the 
radiation of free long waves. If the barred region is sufficiently 
wide, nearly complete reflection of short waves is possible, 
leaving only the free long wave on the transmission side (Hara 
and Mei, 1987). The longshore bars can, therefore, act like a 
magic wall that separates a concert hall: violins played on one 
side sound like cellos on the other. An experimental confirma- 
tion of the theory is shown in Figure.3(from Hara and Mei, 
1987). More research is definitely needed to couple the dynam- 
ics of the sediments with that of the waves. 
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Figure 3. 


Evolution of a group of waves arriving from x < 0 and Bragg- 
reflected by periodic bars extending from x = 0 to 12 mona 
horizontal seabed. The experimental data are shown in solid 
curves; the theoretical wave envelopes are shown by dashed 
curves. Note the dispersion of the incident wave envelope A 
and the reflected wave envelope B. Also shown are the second- 
order long wave trace ¢... From the record a x = 12 m, the 
long wave arrives earlier than the short wave envelope, 
implying the presence of the free long wave. 
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Modulation and instability of short 
gravity waves riding on a long wave 


For weak interactions between two wavetrains of different 
lengths and low amplitudes, Longuet-Higgins and Stewart 
(1960) have shown that short waves become steeper near the 
long-wave crests and flatter near the long-wave troughs. In 
developing their theory, Taylor expansion about the horizontal 
mean surface is made, thus the contrast in amplitudes and 
wavelengths between the two wavetrains must be limited. On 
a finite amplitude long wave, the modulation of a train of 
collinear short waves of infinitesimal amplitude has been 
studied by Longuet-Higgins (1987), who calculated the long 
wave by numerical means. Henyey, et al. (1988) used a Hamil- 
tonian formalism and deduced the action conservation princi- 
ple for infinitesimal short waves. Zhang and Melville (1990) 
went beyond linear modulation to account for higher-order 
effects in short waves and obtained a nonlinear Schr "odinger 
(NLS) equation for the short-wave envelope. Some coeffi- 
cients in the NLS equation have slowly varying coefficients 
representing the effects of long waves. They used the result to 
find the envelope of a Stokes wave, which remains steady with 
respect to the long wave. 

As was pointed out by Longuet-Higgins (1987), an effect 
of the long wave is to provide an oscillating gravity field for 
the short wave. Standing waves in a vertically oscillating basin 
are known to exhibit subharmonic instability. Similar instabil- 
ity must therefore be expected in short gravity waves riding on 
along wave. Naciri and Mei (1992) first re-examined the linear 
effects of a long wave on the modulation of infinitesimal short 
waves over a time span of a few long-wave periods, and then 
studied the sideband instability as well as the subsequent 
nonlinear evolution of Stokes waves over many long-wave 
periods. Since sideband instability is a special case of quartet 
resonance discovered by Phillips (1960), this example gives 
some insight on the effects of a long wave on resonant inter- 
actions. 

Let k, @ and A be the wavenumber, frequency and ampli- 
tude of the short wave and K,Q and B be the corresponding 
quantities of the long waves. Then for infinitely deep water, 
Q=VeK and w=Vgk. If the wave frequency contrast is 
Q/w@ = O (€) << 1, then the wavenumber contrast must be 
K/k= O (e”). Moreover, if the short waves are gentle but the 
long waves are steep, kKA=O(e) << 1 and KB=O(1), then 
kB = O (i/e) >> 1, so that the short wavelength is much less 
than the long-wave amplitude. A Taylor expansion cannot be 
made from the horizontal mean sea level but must be made 
from the long-wave surface. Instead of using orthogonal coor- 
dinates conforming with the long-wave surface, Lagrangian 
coordinates can be used, where the fluid motion is described 
by the instantaneous positions of all fluid particles: 
X(a,b,c,t),Y(a,b,c,t),.Z(a,b,c,t), where (a,b,c) refer to the initial 
position of a particle. The usual kinematic boundary condition 








Figure 4. 


Linear growth rate of sidebands as a function of short wave 
parameter a and the sideband wavenumber v for a long wave 
with KB = 0.3. The short wave steepness is KA = 0.13. 








that a fluid particle on the free surface must never leave the 
free surface can be stated simply by designating the free 
surface as c=0. To avoid excessive computation, Naciri and 
Mei (1992) further modeled the long wave by Gerstner’s 
simple but exact two-dimensional solution 


x=a+B sin (Qt-Ka)eX , 2=c+Bcos(Qt-Ka)eX (2.2) 
with Q?=gK. It is well known that the Gerstner wave has finite 
vorticity that attains the maximum at the free surface c=0 and 
is counterclockwise if the wave propagates to the right; there- 
fore, it cannot be the result of wind blowing in the same 
direction of the wave. 

For a nearly sinusoidal short wave superimposing on 


Gerstner’s long wave, the evolution equation reads, in normal- 
ized form, 


0A 1 
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where 


a= (kKAY* , (2.4) 


pm 
28 Q 


withA , o , k being the characteristic amplitude, frequency and 
wavenumber of the short waves, and 


@=Q1-Ka. (2.6) 
The normalized coordinate @ can be as large as O (e~') , where 
€+0 (Q/w ,kA)<<1. 

For @=0(1), i.e., with a few long-wave periods, all O(e) 
terms can be ignored as a first approximation 


0A __ —KBsing 


= ; 2.7 
a9 1-—2KBcos@ (2.7) 
which can be solved to give 


1Al _ u 
— 2KBcos@ + K*B* 





i (2.8) 





Figure 5. 


Typical evolution of short wave envelope riding on a long wave. 
(a) The Fourier domain where A, is the amplitude of the carrier 
wave, A, is the amplitude of the sideband, and A., A,... are 
the higher harmonics of the side bands. (b) The physical 
domain while gis the time and & the spatial coordinate moving 
at the group velocity. Note the spatial-temporal chaos. 
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Figure 6. 


Dispersion relation for capillary-gravity waves on a current. 
For fixed w and U, the wavenumber x is given by the 
intersection of the straightline w - KU and the curve o(x). In 
Fig. 6.a, there is only one interaction at x,. In Fig. 6.b, there 
are three intersections at x,, K, and K,, with k, < kK, < K, 
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This simple formula agrees well with Longuet-Higgins’ 
(1987) numerical solution of irrotational Stokes wave with 
KB up to 0.3. 

The Stokes waves envelope A , (@) can be regarded as the 
solution of (2.3) independent of € and can be obtained explic- 
itly. Sideband perturbances are governed by a linear partial 
differential system with periodic coefficients and can be stud- 
ied with the aid of the Floquet theory. 

As a sample result from Naciri and Mei (1992), we show 
in Figure 4 for kA = 0.13 ,Q/o =0.07, which corresponds to 
a=0.135, the effects of KB and the sideband wavenumber v. 
When the long wave is absent, KB=0, the band of instability 
reduces to that of Benjamin and Feir (1967) with the band- 
width given by 0 < v <2 ¥2. When KB increases, new bands 
of instability appear and the maximum growth rate increases. 
When the sidebands of the fundamental harmonic are unstable, 
their higher harmonics may also easily fall into one of these 
new bands of instability. Great complexity in the nonlinear 
stage can, therefore, be expected. A sample numerical result in 
Figure 5 (from Naciri and Mei, 1992) shows spatio-temporal 
chaos. As ot increases, the trend toward chaos amplifies. Thus, 
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by the inherent nonlinearity on the ocean surface, long-wave 
modulation of narrow-banded short waves can contribute to 
randomness with scales much greater than the characteristic 
wavelength. This mechanism is, therefore, complementary to 
the randomization of smaller scales due to wind input and/or 
turbulent mixing. 

Instability of collinear short waves on a long wave due to 
oblique sidebands has also been studied by Naciri and Mei 
(1994), but nonlinear evolution has not been examined. A 
quick turn to spatio-temporal chaos is expected, which renders 
the assumption of slow modulation inaccurate. Other theoreti- 
cal approaches allowing broad bands are needed. 


Reflection of capillary-gravity waves 
by a nonuniform current 


When short capillary-gravity waves are reflected by a long 
wave or a nonuniform current, striking physics may arise. In 
a flume with a sloping bottom, Pokazeyeu and Rosenberg 
(1983) and Badulin, Pokazeyeu and Rosenberg (1983) first 
produced a steady and nonuniform current with speeds in the 
range of 4-24 cm/s flowing from shallow to deep water. In the 
direction of intensifying current (from deep to shallow water), 
wave packets of 2-11 Hz carrier frequency were sent. Along 
with complete blockage by the opposing current, the waves are 
found to be shortened by orders of magnitude. Because shorter 
waves are more easily damped by viscosity, the much short- 
ened waves can be completely annihilated without the help of 
breaking (Phillips, 1981). Badulin, et al. (1983) further ob- 
served the phenomenon of double reflection with drastic short- 
ening after each reflection. 

As explained by Badulin, et al. (1983), the shortening can 
be understood by examining the dispersion relation. For an 
opposing current with local speed U < 0, the dispersion rela- 
tion is 


"% 
o+kiui=(gk+ 2 =o(k). 
\ p 


In Figure 6, the function o(k) is represented by the curve, while 
the left-hand side by the straight line with the slope 
@/k=\U\. Their intersection gives the wavenumber of the 
short wave. Let @, be the intersection of the straight line 
tangent to the o curve at the point of inflection. We can 
distinguish two cases. 

(i) @ > @,: As shown in Figure 6.a there is only one 
intersection at the wavenumber k, >0 which increases with 
\U\. Hence, there is only one wave propagating into the current. 
As it propagates forward to the shallow water, waves becomes 
shorter. Since do/dk >—U ,C, +U>0,-U,C,+U 0, energy 
is propagated forward into shallow water also. 





(ii) @< @,; In the region where IUI < IU,I, there is only 
one solution k,. Again, do/dk + U = C,+U > 0, signifying for- 
ward propagation of wave phase and energy in the current. 

As the wave enters the zone of stronger current 
IU,! < IU <1U,I, three waves coexist: 0 < k, < k, < k, (Figure 
6.b). For k, and k, do/dk >— U, so that the energy is trans- 
ported forward. For k,, however, do/dk > — U, hence energy 
is reflected. As |U\ further increases beyond 1U)], i.c., 
U1 >1U,I, only the shorter wave (k,) can exist with energy 
transported forward. 

In summary, there are two points of reflection: one at x, 
(where U=U,) and one at x, (where U=U,,). Between them k, 
can be regarded as the incident wave, k, the reflected wave, 
and k, the transmitted wave. The wavelength is greatly reduced 
after each reflection. 

This problem can be analyzed by the uniformly valid 
theory of Shyu and Phillips (1992). Trulsen and Mei (1993) 
employed a WKB approximation away from, and boundary- 
layer corrections near, the reflection points. With a crude 
account for viscous damping, the dramatic shortening and 
effective annihilation of a long-wave packet after double re- 
flection can be clearly seen in Figure 7 (from Trulsen and Mei, 
1993). 


Triad resonance of capillary-gravity 
waves affected by a long wave 


At the beginning of this century Harrison (1909) and 
Wilton (1915) found that a monochromatic progressive capil- 
lary-gravity wave can excite its second harmonic to significant 
amplitude. Half a century later McGoldrick (1965) found this 
phenomenon to be a special case of nonlinear triad resonance, 
where three capillary-gravity wavetrains can resonantly ex- 
change energy among one another through quadratic coupling. 
After the seminal works by Phillips (1960) on gravity waves 
and by McGoldrick (1965) for capillary-gravity waves, the 
study of nonlinear resonance has mushroomed. For a compre- 
hensive survey, see Craik (1985). 

The effect of long waves on a resonant triad of cap- 
illary-gravity waves has been treated recently by Trulsen 
and Mei (1995), who also used Lagrangian coordinates 
and made the following order assumptions: 
kA = O() KB = O(€) K/k = O(€”)Q/w = O(e). With these 
scale choices, both short and long waves are gentle in slope. 
By ignoring viscous damping and applying the multiple-scale 
analysis to the second order, three evolution equations are 
derived: 


dA, ; — 
> +C,1- VA; + iB,Ajcos@ + io1,A7A3 + 0 


y (2.10) 


(2.11) 


dA, St 
>= C,; - VA2 + i,Agcos@ + iaAjA; +0 


dA; . 7 
> +s - VA; + iB,A;Cosp + i0,AjA,+0 (2.12) 


The terms representing interactions between short and long 
waves are B; A,cos@, where B; is proportional to the long-wave 
amplitude and @=—Q is the long-wave phase. Trulsen and Mei 
(1995) focused their attention to the case with spatial uniform- 


ity so that an ordinary dynamical system of the sixth order is 
obtained. By a judicious transformation and the use of conser- 





Figure 7. 


Snapshots of waves advancing into an opposing current. From 
the top figure to the bottom, time increases. Damping is 
included. Note the annihilation of waves after second reflec- 
tion. 




















































































































Twol1996 9 





vation laws, the number of independent variables is reduced to 
two, yielding 


=-r{(l1-pj(i+j,)] 2 sing (2.13) 


and 


dy 37 -2j+ji-i. 
— =rcosy e4s » \71.2 
dt 211 -pj G+.) 





+A+6 cos t (2.14) 


where j is proportional to the wave action of the shortest wave 
in the triad, and y=, +@,—@, is a special combination of phase 
angles of the wave actions. Because of the omission of viscos- 
ity, the system is Hamiltonian. Without long waves the above 
Hamiltonian system is autonomous and may be solved exactly. 
For weak long waves, a perturbation analysis is possible, 
which shows various possibilities of bifurcations and modula- 
tional resonances of superharmonic or subharmonic varieties 
at certain long-wave frequency. As the long-wave amplitude 
increases, these resonances lead to homoclinic orbits, horse- 
shoe tangles, and chaos, as is expected of a Hamiltonian 
system. A sample result is shown in Figure 8 (from Trulsen 
and Mei, 1995). This study shows again that a long wave 
triggers chaos among nonlinearly resonating waves. Further 


studies of equations (2.10-12) on the spatio-temporal chaos are 
worthwhile. 


Frequency downshift and 
long-time modulation 


The prediction of the spectral development of wind-in- 
duced waves is the primary goal of wave-modeling, which is 
immensely important to engineering and oceanography. In this 
task one must integrate the mechanisms of wind input, 
wave/wave interactions and dissipation due to breaking, inter- 
nal and bottom friction, etc. Based on existing theoretical, 
laboratory and field knowiedge of these mechanisms, impres- 
sive progress has been made, resulting in the third-generation 
model for wave prediction (Koman, et al., 1995). On the other 
hand, a deeper understanding is also being pursued by studying 
separately various aspects of the whole phenomenon. 

One prominent feature of a wind-wave spectrum evolu- 
tion is called the frequency downshift. Under stationary off- 
shore wind, the frequency at the peak of the energy spectrum 
of wind-generated waves is known to steadily decrease with 
increasing fetch, while the peak height increases. From JON- 
SWAP measurements, the peak frequency varies with fetch 
according to 


£. = 3.5 (gx/U 49) ~033 (3.1) 
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where x denotes the fetch, and U,, denotes the mean wind 


speed 10 m above the mean sea level. The cause of the 
downshift of peak frequency has been the focus of attention 
by various authors in recent years. 

In their experimental study of the Benjamin-Feir instabil- 
ity of uniform Stokes waves, Lake, et al. (1977) sent steep 
waves down a flume of 40-foot length. For Stokes waves of 
relatively low initial steepness, unstable sidebands grew at 
equal rates. The time scale of growth decreases with the square 
of wave steepness. However, for steeper waves the most un- 
stable sidebands grew at equal rates only near the wavemaker. 
Further downstream the lower sideband becomes higher than 
both the upper sideband and the carrier wave. As a conse- 
quence, the spectral peak shifts to a lower frequency. Without 
breaking, there is a tendency for the carrier to regain its 
strength and for the sidebands to diminish further down the 
tank. Melville (1982) performed similar experiments in a 
longer flume with steeper waves. He observed not only faster 
growth of the lower sideband but wave breaking at the height 
of modulation. 

For low enough steepness, the first nonlinear stage of 
sideband evolution is well described by the nonlinear Schr 
"odinger equation (NLS) 


0A a°A 


2 
= p Low 4 Bes 
is, tes5 +B ul 0 (3.2) 


Because of the even derivatives with respect to x, the nonlinear 
Schrodinger equation predicts that the sidebands, if initially 
equal, will grow and decay at equal rates. Using a numerical 
spectral method to solve the integral equation of Zakharov, 
which has a better accuracy than NLS, Yuen and Lake (1980) 
showed that the lower sideband grows faster than the upper 
sideband, during the initial phase of evolution. After the mo- 
ment of maximum modulation, the evolution depends on the 
initial wave steepness. For initially gentle waves, energy in the 
sidebands is returned to the carrier wave, and the process of 
modulation and demodulation remains cyclic. For steep 
waves, the evolution can become chaotic. These numerical 
results appear to have been computed for wave-steepness too 
large to be consistent with the theory, however. 

The instability of infinitesimal sidebands has also been 
studied for a steep wave by Longuet-Higgins (1987). By 
adding fourth-order terms to the nonlinear Schr "odinger equa- 
tion, Dysthe (1979) has brought the predicted rate for an 
unstable growth rate closer to the more exact numerical results 
of Longuet-Higgins. The extension of NLS equations reads: 


a 2kox = 8k ax 


+5w \Al7A 








Figure 8. 


Typical Poincare maps of the resonant triad under the influence 
of long wave. (a) B = 0.1, (b) B= 0.2 
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where @ represents the potential for a wave-induced current. 
Because of the presence of odd derivatives at the fourth order, the 
growth of unstable sidebands should no longer be symmetric after 
a long enough time of the order k* a*ax=0(1). By a numerical 


method Lo and Mei (1985) solved (3.3) to (3.5) for a Stokes 
wave with a small periodic modulation corresponding to the 
most unstable sidebands. They found that the initial stage of 
instability is followed by the faster growth of the lower side- 
band, which may overwhelm the carrier waves at the height of 
modulation. After some time energy returns from the side- 
bands to the carrier wave. This modulation-demodulation 
process goes on periodically over the time scale of 
«t(ka)"=O(1). Lo and Mei (1985) also introduced viscous dissi- 
pation damping in the boundary layers along flume walls. The 
recurring trend is qualitatively unchanged. Applications of these 
equations are more successful in the prediction of the disintegra- 
tion of an isolated wave packet into a string of smaller packets as 
seen in the experiments by Feir (1967). 

In the flume experiments by Melville (1982), breaking is 
the most pronounced at the time when the sidebands attain 
their maximum. To simulate this effect, Trulsen and Dysthe 
(1990) introduced a simple model simulating energy loss by 
breaking, as well as wind input (1992). These two effects are 
represented by a pair of source and sink terms in (3.3) to (3.5). 
The source term is proportional to the product of wind stress 
and the wave amplitude. The sink term is of the relaxation type, 
which ensures energy loss when the wave amplitude exceeds 
a certain ceiling. At zero wind speed, breaking is found to 
preserve the dominance of the lower sideband. At higher wind 
speed, the higher harmonics of the lower sideband may main- 
tain further downshift. For a still stronger wind, the carrier 
waves attain a critical level beyond which sidebands do not 
grow at all; they are overdamped. While qualitatively encour- 
aging, the model depends on coefficients that can only be 
determined by curve-fitting with measurements yet unavail- 
able. Recently, Okamura (1995) studied the evolution of stand- 
ing waves and introduced two models for breaking. He notes 
from the experimental records of Melville and Rapp (1988) 
that at breaking the decreases in local velocity are much more 
pronounced than those in the local free-surface height. In his 
two models, the kinetic energy is removed when the sidebands 
are at their peak. The rate of removal is represented by reduc- 
tion factors that are not deducible by theoretical arguments. 

In amore deductive theory, Hara and Mei (1991) analyzed 
the downshift problem by considering weak wind and internal 
dissipation, but no breaking. The basic wind profile is linear 
but the wave-induced disturbances in air are treated as pseudo- 
laminar (constant eddy viscosity). By analyzing the wave-in- 
duced boundary layers in both air and water, they found that 
the Doppler shift caused by the wind-induced drift near the 
water surface must be important. The predicted effects of wind 
on sideband instability is reasonably well supported by the 
flume experiments of Bliven, et al. (1986). Furchermore, a 
monotonic trend of frequency downshift is found in the non- 
linear stage (Figure 9, from Hara and Mei, 1991). Thus the 
combined effect of wind and internal dissipation can be a 
possible mechanism to maintain the frequency downshift. 
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In a subsequent paper Hara and Mei (1994) studied the 
linear and nonlinear phases of sideband instability of capillary- 
gravity waves under wind and dissipation. They found that for 
waves longer than 27/k , there is a frequency downshift after 
a long time. But for waves shorter than 27/k , there is an 
upshift instead. The critical number corresponds k, to the 
capillary-gravity wave with the lowest phase speed and to the 
resonance of Wilton’s ripples. 

For sufficiently strong winds, internal friction is too weak 
to counterbalance the trend for growth; waves must break. A 
consistent theory must account for strong nonlinearity and 
turbulence in breaking waves; this task will remain a challenge 
to the wave scientist for some time. 

We now turn to some recent studies that have strong 
engincering implications. 


Nonlinear wave problems of 
engineering interest 


Envelope soliton in a ship wake 


Very long and narrow ship wakes of V shape have been 
revealed as bright lines in satellite photographs (Fu and Holt, 
1982; Shemdin, 1990). The half angle of these wakes are 
typically 2°-3°, hence much narrower than the classical Kelvin 
wake with a half angle of 19.7°. The length of such wakes can 
be 15-20 km. Since the relation between ships and their wakes 
is of obvious interest to aerial surveillance, numerous expla- 
nations have been proposed for these V wakes. For example, a 
point disturbance moving steadily in a two-layered ocean is 
known to generate interfacial Kelvin waves with an apex angle 
smaller than 19.7°. It therefore has been suggested that these 
Kelvin waves can modulate the short gravity-capillary waves 
on the sea surface and give rise to the observed V-wake pattern 
(Tulin and Miloh, 1990). However, the modulated surface 
wave may be too weak to cause such a visible radar backscatter 
for a mild wind and a mixed-layer depth typical in the ocean 
(Shemdin, 1990). Another possibility is the incoherent oscil- 
lations associated with turbulence in the wake (Munk, et al., 
1987). Yet another explanation attributes the V wake to the 
surface manifestation of vortices or bubbles generated from 
the ship’s propellers. A recent theory by Gu and Phillips (1995) 
reasoned that the coherent structure in the turbulent jet gener- 
ated by the propellers may produce oscillations along the jet 
boundaries, which act like oscillating pistons along the vertical 
center plane of the wake. These oscillations radiate two groups 
of small-scale surface waves that are seen as the two arms of 
the V wake. 

The dominant mechanism of the V wake may be different 
in different observations; a single explanation may not fit all 
observations. To provide more definitive data, Brown, et al. 
(1989) conducted field tests near the coast of Southern Cali- 
fornia. Along the path of a 25-m coastguard cutter cruising at 
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the speed of 7.7 m/s, detailed records were taken at stationary 
wave masts. Instead of the very narrow V wake, their records 
showed strong evidence of envelope solitons consisting of 
waves of frequency centered around 3.28 rad/sec. These enve- 
lopes formed V-like arms with a half angle 11 degrees from 
the ship’s path, and extended roughly 1 km behind the ship. It 
is unclear whether these soliton wakes correspond to the much 
narrower and longer V wakes that motivated the field test; they 
nevertheless warrant explanation in their own right. 

Much earlier, Lighthill (1967) predicted by a linear theory 
that a point pressure traveling on the free surface and oscillat- 
ing at a constant frequency generates a Kelvin wake with a 
smaller apex angle than a steady Kelvin wake. Suggested by 
the comments of Brown, et al. (1989) that wind in the speed 
range of 0-6 m/sec and wind-induced swell of 0-3 m height 
might have been present during their field tests, Mei and Naciri 





Figure 9. 


Evolution of carrier wave amplitude and the sidebands under 
the influence of wind and dissipation. In Fig. 9.a, B, is the carrier 
wave. B , the lower sideband, and B,, the upper sideband. In 
Fig. 9.b, B,, are the second harmonics of B.,. Fig. 9.c gives the 
spectrum of wave envelope at dimensionless time t = 50. 












































(1991) constructed a theory for a slender ship that oscillates 
while cruising. The ship acts as an oscillating piston advancing 
at a constant speed, and radiates short plane waves away from 
the ship. The linear part of the mechanism is therefore partially 
similar to that proposed by Gu and Phillips (1995). Viewed 
from a coordinate system fixed on the ship, the incoming head 
sea propagates in a current of speed V. By Doppler’s effect, the 
ship oscillates at the encounter frequency given by 


Ve 
gcosa |” 


where @, is the frequency of the ambient waves in a fixed 
coordinate system. For a head sea o=r,, 


Vo, 
@ = Wo —S >@,; 


The records of Brown, et al. (1989) suggest that a head sea of 
@, = 1.5 rad/sec would give rise to a head sea of @ = 3.28 
rad/sec. The ship can be induced to oscillate at the encounter 
frequency. 

Oblique soliton envelopes are known to be solutions of 
the two-dimensional nonlinear Schrodinger equation. Rewrit- 
ten in a moving coordinate fixed on the ship, the envelope 


equation is 
vile ta od + « ll? A=0, (4.3) 


where the ship’s path coincides with the y-axis and current 
speed V is in the positive y direction. Let e=kA,<<1. We assume 


the ship to be slender so that kKB=O(1) and kL= O (1/€)>>1 
where B and L are the beam and length of the ship. Very close 
to the ship kx< O(/e) ky=O (1), waves generated by the 
oscillating and advancing ship are described by the linear 
theory. In the intermediate domain kxky = O(1/e), the evanes- 
cent modes have lost their significance; only propagating 
waves remain. These outgoing waves are nearly plane waves 
with a finite crest length which is much greater than the 
wavelength, hence are slowly modulated in y and obey the 
linear envelope equation given by the first three terms in (4.3). 
In the far field where kx= O (1/e”) , ky = O(1/e), nonlinear 
effects become important and we must use (4.3) with the 
boundary condition that A(O,y) is given along x=0. Trans- 
formed to a moving coordinate, the nonlinear Schr "odinger 
equation can be reduced and normalized to 


-ite z[1-(F [Ss +era-o, 


(4.1) 


(4.2) 


Ho MM, 
ot ax = y 


(4.4) 


where a= kA, &=K A, (kx - (Cg/V) ky ),and 1) =k A, (Co/V) ky. 
By a known property of the nonlinear Schrodinger equation, 
solitons exist only if the coefficients in the square brackets 
above are positive, 


V Cy (4.5) 


The apex angle of the classical Kelvin wave is defined by C 
g/V <1/V8 . Thus, if 


1 Gc 1 
8 Vv <v2 


solitons exist and lie outside Kelvin’s wake. Otherwise, if 


(4.6) 


Ce 1 
V < Te (4.7) 


solitons lie inside Kelvin’s wake. Since for deep water waves 


(4.8) 


a small C,/V or a high-frequency oscillation gives rise to a 
small wake angle. 

Mei and Naciri (1991) simulated the effect of heaving and 
pitching by adopting appropriate profiles for A as an initial 
condition for (4.4), which is then solved numerically. After a 
transfomration of coordinates from & to o where 


53h 


“Fane (4.9) 


with &) being proportional to the ship length L and the wave 
steepness KA , 


= RAL 


(4.10) 


Equation (4.4) becomes 


j24,1 70, apg no 


-is +359 (4.11) 


For a heaving ship the inital envelope close to the ship is 
modeled qualitatively by 


a(S ,0)=sech(G/6p) (4.12) 


This profile represents a solition if 6,=1/V2. . If the ship is 
longer or if the oscillation amplitude is higher,the value of 6, 
is greater. In Figure 10 the spatial evolution of the initial 
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Figure 10. 


Evaluation of soliton envelope in the wake of a heaving ship 


for various values of 60; oo = 1.6 (Fig. 10a), 3 (Fig. 10.b), and 
5 (Fig. 10.c) 
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envelope on one side of the ship path is shown for several 65. 
More solitons appear for larger 6. The number of solitons is 


reflected in the number of undulations in the envelope. 
Whether these undulations correspond to the observations is 
not easy to decide on the basis of the available data. It would 


be interesting to check these predictions in controlled experi- 
ments in a large basin. 


Upstream radiation of solitons 
in a channel 


Nonlinear shallow water waves have been the subject of 
intensive scientific study in the past three decades. The seminal 
survey by Ursell (1953) on the long-wave paradox resolved 
the differences between the theories of Airy and of Rayleigh 
and Korteweg-de Vries. In the former, any one-dimensional 
compressional wave must break because the governing equa- 
tions are hyperbolic. In the latter, smooth waves of permanent 
form can exist. The subsequent flood of activities owes much 
to the numerical experiments of Zabusky and Kruskal (1965) 
and to the use of the inverse-scattering theory to solve an 
initial-value problem for the Korteweg-de Vries (KdV) equa- 
tion exactly. These developments have been extensively re- 
viewed in numerous articles and books, e.g., Ablowitz and 
Segur (1981). 

A two-dimensional extension of the KdV equation is the 
K-P (Kadamtsev-Petviashvilli) equation 


(u, + Guu, + Ux), — uy, =0, (4.13) 
which holds that if the characteristic lengths in x and y direc- 
tions are related by O(L,) = O(L2), numerical and analytical 
solutions of (4.9) have revealed the mechanics of Mach stems 
(Miles, 1977; Chen and Liu, 1988) and the existence of two- 
dimensional periodic waves (Hammack, et al., 1989). 

An unusual behavior of ships in restricted water was noted 
in model tests before WWII by Thews and Landweber (1935) 
in the U.S. and Izubuchi and Nagasawa (1937) in Japan.* 
When a model ship is towed near the critical speed, not only 
does the ship residence suddenly become very large but the 
flow around the ship is unsteady. Here is a vivid description 
by Kinoshita (1946): 

"When the ship speed is higher than V’ the whole quantity 
of water flowing against it can no longer pass, but .... is 
dammed up in front of the ship and consequently turned into 
backwater moving ahead of it." 

Later, Graff, et al. (1964) had a similar experience in their 
experiments, and reported the measured ship resistance only 
after averaging over the time-dependent records. 

Before the age of nonlinear instabilities and resonances, 
steady inputs were thought to produce only steady responses. 
Thus several authors extended the steady transonic aerody- 
namic theory to ship waves (Maruo, 1948; Lea and Feldman, 
1972; Maruo and Tachibana, 1981). It was finally realized 
experimentally by Huang, et al. (1982) and numerically by Wu 
and Wu (1982) that the waves generated by a ship advancing 


* Research in Japan before WWII was prompted by ship operations in the rivers of Chinese Manchuria under Japanese occupation. (Private communication 


by Professor H. Maruo of the Yokohama National University, 1985). 
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near the critical speed are inherently unsteady. Unlike gas 
dynamics, water waves are distinguished by their dispersive- 
ness. Thus by accounting for unsteadiness, dispersion and 
nonlinearity, theories for two-dimensional disturbances were 
advanced by Wu and Wu, 1982; Akylas, 1984; Ertekin, et al., 
1986; and Wu, 1987. The phenomenon is governed by an 
inhomogeneous KdV equation. 

Replacing the near field of a ship hull by an equivalent 
wall-sided strut of width B, Mei and Choi (1987) showed the 
equation for the far field of a slender ship to be 


2-2m 
ve 3€ 1 
Pup =— 7 |e 2-32 +3 Pa » @.14) 


where all quantities are dimensionless with x normalized by 
the ship length L and y by the channel width W so that h=“v. 
The dimensionless ratios are, respectively, 


L 
€=A/h ° w=h/L, 10 

and @ is the dimensionaless measure of detuning from the 
critical speed. The boundary conditions are 


1 


B 
%,=7 Y, on 7 =0 
Ley * 


(4.15) 


along the hull and 


9, =9 on H=1 (4.16) 


along the channel bank. The channel width is assumed to be 
W/L=O(u™) so that )=O(1). Thus, for a wide channel so 


that m=1, (4.10) is just the K-P equation, which can be solved 
numerically. For a narrow channel with m= 1/2, say, Pan 
vanishes to the leading order, thus the wave field is approxi- 
mately one-dimensional fore and aft. Eqn. (4.10) can be inte- 
grated to give an inhomogeneous KdV equation 
3 1 1 
Ge ~ 08 — 5 SE. — § Sax = BY (4.17) 
where E=-9.. 

Mei (1986) used (4.13) to predict the one-dimensional 
upstream solitons measured by Ertekin and found good agree- 
ment despite the blunt bow of the model ship and the two-di- 
mensional waves in the wake. By using matched asymptotics, 
Mei and Choi (1987) deduced the equivalent strut width for a 
ship hull with a finite draft, and calculated the oscillatory 
forces on the ship due to the radiation of solitons. Recently, 
Chen and Sharma (1995) have improved the numerical scheme 
for solving the K-P equation and the accuracy of the ship hull 
boundary conditions. A sample result of their computation is 


shown in Figure 10 (from Chen and Sharma, 1995), showing 
the one-dimensional solitons radiated upstream and the two- 
dimensional waves downstream. 

Besides engineering interests, the problem of upstream 
solitons has been studied in the oceanographic context. For 
example, interfacial waves in a two-layered sea can be induced 
by tides passing a topography. Nonlinear solitons can be 
expected to appear in a strait if a long topographical feature is 





Figure 11. 


Free surface evolution due to a ship advancing at the critical 
speed in a shallow channel. As time passes, more solitons are 
radiated upstream. (Chen and Sharma, 1995). 
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present either across the strait or along a coast (e.g., Melville 
and Helfrich, 1987). In a two-layered ocean two linear modes 
of internal waves with different phase speeds can exist. When 
the tidal flow speed is near either wave speed, upstream 
solitons are possible. For a continuously stratified sea, there 
can be many internal wave modes. Nonlinear dispersive effects 
near one of the critical wave speeds can be similarly complex 
and deserve further study. 


Subharmonic resonance of 
articulated structures 


The fabled city of Venice is situated in a small lagoon on 
the shore of the Adriatic Sea. Due to the dense population and 
industrialization, excessive use of groundwater has led to the 
lowering of land surface. Storm surges, which enter the lagoon 
through three inlets, now flood the city as often as 50 times per 
year. To save the ancient city from deterioration, engineers in 
Italy have proposed to construct a storm barrier across each of 
the three inlets (Bandarin, 1991). To allow tidal flushing of the 
lagoon and to preserve the scenic beauty, the barriers must be 
kept out of sight except during a storm. A design has been 
under consideration since 1984 for a barrier consisting of many 
steel boxes hinged on acommon axis that is fixed on the seabed 
across each inlet. These boxes are unconnected with each other 
in order to reduce the cost for a strong support and for easy 
maintenance. When a storm is forthcoming, air is compressed 
into the boxes so that they rise to an inclination of 60 degrees 
from the horizontal. Because the inlet channels are long, the 
storm waves are expected to arrive with their crests parallel to 
the gates. Then the gates are expected to swing to and fro in 
unison, transmitting relatively little force and torque to the axis 
while serving as a dam to stop the ocean water from entering 
the lagoon. 

Laboratory tests in Delft and in Italy have revealed an 
unexpected phenomenon, however. When the incident waves 
arrive at certain frequencies, the gates do not swing in unison. 
Instead, neighboring gates oscillate out of phase, therefore 
diminish the intended function as a dam. These laboratory tests 
also showed that the out-of-phase oscillations have frequen- 
cies just one-half of that of the incident waves. 

Decades ago, subharmonic resonance was studied ac- 
tively for an unrelated oceanographic phenomenon - surf beats 
that may contribute to the formation of beach cusps. They can 
also cause hazards in resonating harbors or for moored ships. 
It is known that along a sloping beach there can be natural 
modes called edge waves with energy trapped along the shore. 
Guza and Bowen (1976) showed that an edge wave can be 
excited subharmonically by normally incident wave. Theoreti- 
cal analyses for steady incident waves were advanced by 
Rockliff (1978) and Miles (1990). 

Mei, et al. (1994) found that the articulated construction 
of the Venice gates can support trapped waves that do not 
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radiate energy to infinity. In fact, there can be a discrete set of 
modes where neighboring gates oscillate out of phase without 
radiating wave energy. They have tested the existence of one 
trapped mode by comparing the eigenfrequencies with experi- 
ments for two vertical gates in a narrow wave flume. Moreover, 
by sending normally incident waves at twice the eigenfre- 
quency of a trapped mode (see Figure 11), resonance was 
produced. P. Sammarco has recently confirmed that the non- 
linear behavior of the gate displacement B is governed by the 
Landau equation 


~iF = 85+ oAB* +B BPE, 


(4.18) 


where A is the incident wave amplitude. Mei and his colleagues 
P. Sammarco and H. Tran are now examining the stability, 
bifurcations and possible chaotic motions for periodically modu- 
lated incident waves with A = A, (1 + 1 cosQ2), where Q is a 
very small modulational frequency. An advance in modern sci- 
ence is finding its usage in a traditional engineering problem that 
dates back to 4000 BC, when Egyptians built the first breakwaters 
for the harbor of Alexandria. 


Concluding Remarks 


From the physical point of view, one of the most striking 
lessons learned in recent advances is that, given long enough 
time or large enough space, even weak nonlinearity can lead 
to significant changes to the linearized theories. Thus steady 
forcing can result in a transient response; deterministic input 
can yield a chaotic evolution. Many long-cherished concepts 
of linearized theories have fallen by the wayside. While non- 
linearity is drawing attention in many scientific disciplines, the 
dynamics of ocean surface waves is still the most furtile 
ground for studying its impacts and finding more challenges. 

From the viewpoint of theoretical analysis, efforts to gain 
new understandings have benefited significantly by the en- 
hanced computer power in numerical simulation. It is, how- 
ever, safe to predict that in the forseeable future analytical and 
numerical approaches will remain complementary for a long 
time; one cannot replace the other. The analytical approach is 
more effective when there are many different scales, especially 
for small amplitudes or weak nonlinearity. Numerical ap- 
proaches are indisputably more powerful when dealing with 
finite amplitudes, but their power diminishes when many 
different scales and two or three dimensions must be consid- 
ered. Future progress is likely dependent on the combination 
of the two approaches in some optimal way. 


Acknowledgment 


The author acknowledges with gratitude the sustained 
ONR support received since the beginning of his academic 





career (Program of Fluid Mechanics, Program of Ocean Engi- 
neering, and The Accelerated Research Initiative on Nonlinear 
Ocean Waves). Partial support from the National Science 
Foundation for many years is also recorded with pleasure. 


Biography 


Chiang C. Mei is Edmund K. Turner Professor of Civil 
and Environmental Engineering at Massachusetts Institute of 
Technology. He received his PhD from Caltech in Engineering 
Science. After a two-year post-doctoral fellowship there, he 
moved to MIT, where he has remained ever since, except for a 
few sabbatical leaves. His research interests are: (i) ocean 
surface waves, including the fluid dynamical theory and appli- 
cations in coastal and offshore engineering, (ii) dispersion of 
suspended particles near, and transport of cohesive fluid mud 
on, the seabed, ‘iii) dynamics of poroelastic media with appli- 
cations to seabed mechanics and land subsidence, and (iv) 
diffusion and transport of volatile organic chemicals in porous 
media. 


References 


1. Ablowitz, M.J. & SegurR, H. (1981). Solitons and the 
Inverse Scattering Transform, Society Industrial and Ap- 
plied Mathematics, Philadelphia. 


Akylas, T.B. (1984). J. Fluid Mech. 141: 455-466. Agnon, 
Y. & Mei, C.C. (1988). J. Fluid Mech. 195: 201-222. 
3. 
Agmon, Y. & Mei, C. C. (1988) J. Fluid Mech. 195: 
201-222 
Badulin, S.I., Pokazeyeu, K.V., & Rosenberg, A.D. 
(1983). Izv. Atmos. Ocean. Phys. 19 (10): 782-787. 
Bandarin, F. 1991 Coastal Zone 91, Amer. Soc. Civ. 
Engrs. 3352-3370. 
Battjes, J.A. (1988). Annu. Rev. Fluid Mech. 20: 257-293. 
Benjamin, T.B. & Feir, J.E. (1967). J. Fluid Mech. 27: 
417-430. 


Bliven, L.F., Huang, N.E., & Long, S.R. (1986) J. Fluid 
Mech. 62: 237-260. 


Bowen, A.J. (1969). J. Marine Res. 27: 206-214. 


. Brown, E.D., Buchsbaum, S.B., Hall, R.E., Penhune, J.P., 
Schmitt, K.F., Watson, K.M., & Wyatt, D.C. (1989). J. 
Fluid Mech. 204: 263-293. 


. Carter, T.G., Liu, PL.F., & Mei, C.C. (1973). J. Water- 
ways, Harbours Coastal Eng. Div. ASCE 99: 165-184 


. Chen, X. & Liu, Y. (1988). Acta Mechanica Sinica, 4 (3): 
201-210. 


. Chen, X. & Sharma, S.D. (1995): J. Fluid Mech. 291: 
263-285. 


. Craik, A.D.D. (1985). Wave Interactions and Fluid Flows. 


Cambridge University Press. 


. Davies, A.G. & Heathershaw, A.D. (1984). J. Fluid Mech. 


144: 419-443. 


. Debnath, L. (1994). Nonlinear Water Waves. Academic 


Press. 


. Dysthe, K.B. (1979). Proc. R. Soc. Lond. 369: 105-114. 
. Ertekin, R.C., Webster, W. C., & Wehausen, J. V. (1986) J 


Fluid Mech. 1659: 275-292. 


. Feir, J.E. (1967). Proc. R. Soc. London A 299: 54-58. 
. Fu, L.L. & Holt, B. (1982). Jet Propulsion Lab. Publ. 


81-120. 


. Graff, W., Kracht, A., & Wwinblum, G. (1964). Trans. 


Soc. Nav. Archit. Mar. Engnrs. 72: 374-401. Gu, D. F. & 
Phillips, O.M. (1995). J. Fluid Mech. 275, 301-321. 


. Guza, R.T. & Bowen, A.J. (1976). J. Marine Res. 34: 


269-293. 


. Hammack, J., Scheffner, N., & Segur, H. (1989). J. Fluid 


Mech. 209, 567-589. 


. Hara, T. & Mei, C.C. (1987). . Fluid Mech. 178: 221-241. 
. Hara, T. & Mei,C.C. (1991). J. Fluid Mech. 230: 429-477. 
. Hara, T. & Mei, C.C. (1994). J. Fluid Mech. 267: 221-250. 
. Harrison, W.J. (1909). Proc. Lond. Math. Soc. 7: 107-121. 
. Heathershaw, A.D. (1982). Nature 296: 343-345. 

. Henyey, F.S., Creamer, D.B., Dysthe, K.B., Schulz, R.L., 


& Wright, J.A. (1988). J. Fluid Mech. 189: 443-462. 


. Huang, D.B., Sibul, O.J., Webster, W.C., Wehausen, J.V., 


Wu, D.M., & Wu, T.Y. (1982). Proc Conf. on Behavior of 
Ships in Restricted Waters, Varna, Bulgaria, 26.1-26.9. 


. Izubuchi, T. & Nagasawa, S. (1937). Jap. Soc. Nav. Archit. 


61, 165-206. 


. Kinoshita, M. (1946). Jap. Soc. Nav. Archit. 76 174-213. 
. Koman, G.J., Cavaleri, L., Donelan, M., Hasselmann, K., 


Hasselmann, S., & Janssen, P.A.E. (1995)/ Dynamics and 
Modelling of Ocean Waves. Cambridge University Press. 


. Lake, B.M., Yuen, H.C., Rungaldier, H., & Ferguson, Jr., 


I.N.E. (1977). J. Fluid Mech. 83: 49-74. 


. Lea,G.K. & Feldman, J.P. (1972). Proc. 8th Suymp. Naval 


Hydrod Paris, 1527-1651. 


. Lighthill, M.J. (1967) J Fluid Mech. 27: 725-752. 
. Liu, PL.-F. (1989). J. Fluid Mech. 205: 163-170. 
. Liu, P. L-F, (1994) Advances in Coastal and Ocean Engi- 


neering, 1, World Scientific. 125-158. 


. Lo, E. & Mei, C.C. (1985). J. Fluid Mech. 150: 395-416. 
. Longuet-Higgins, M.S. (1970a,b). J. Geophys. Res. 75: 


778-6789; 6790-6801. 


. Longuet-Higgins, M.S. (1987). J. Fluid Mech. 177: 293- 


306. 


Twol/1996 17 





. Longuet-Higgins, M.S. & Stewart, R.W. (1960). J. Fluid 
Mech. 8: 565-583. 


. Maruo, H. (1948). Jap. Soc. Nav. Archit. 84 (1952): 1-15. 


. Maruo, H. & Tachibana, T. (1981). Jap. Soc. Nav. Archit. 
150: 181-187. 


. McGoldrick, L.F. (1965). J. Fluid Mech. 21: 305-331. 
. Mei, C.C. (1985). J. Fluid Mech. 152: 315-335. 
. Mei, C.C. (1986). J. Fluid Mech. 162: 53-67. 


. Mei, C.C. (1989). The Applied Dynamics of Ocean Sur- 
face Wave. World Scientific. 


. Mei, C.C. & Agnon, Y. (1989). J. Fluid Mech. 208: 595- 
608. 


. Mei, C.C. & Benmoussa, C. (1984). J. Fluid Mech. 139: 
219-235. 


. Mei,C.C. & Choi, H.S. (1987). J. Fluid Mech. 179: 59-76. 


. Mei, C.C. & Liu, P.L.-F (1993). Annu. Rev. Fluid Mech. 
25: 215-240. 


. Mei, C.C. & Naciri, M, (1991). Proc. R. Soc. Lond. A 432: 
535-546. 


. Mei, C.C., Sammarco, P., Chan, E.S., & Procaccini, C. 
(1994). Proc R. Soc. Lond. 444: 257-265. 


. Melville, W.K. (1982). J. Fluid Mech. 115: 165-185. 


. Melville, W.K. & Helfrich, K.R. (1987). J. Fluid Mech. 
178: 31-52. 


. Melville, W.K. & Rapp, R.J. (1988). J. Fluiud Mech. 189: 
1-22. 


. Miles, J.W. (1977). J. Fluid Mech. 79: 171-179. 


. Miles, J.W. (1990). J. Fluid Mech. 214: 43-57. Molin, B. 
(1982). Pap. Ship Res. Inst., Tokyo, Japan 68. 


. Munk, W.H., Scully-Power, P. & Zacheriasen, F. (1987). 
Proc. R. Soc. Lond. A. 412: 231-254. 


. Naciri, M. & Mei, C.C. (1992). J. Fluid Mech. 235: 
415-452. 


. Naciri, M. & Mei, C.C. (1994). Wave Motion. 20: 211- 
232. 


. Okamura, M. (1995). Wave Motion. 
. Phillips, O.M. (1960). J. Fluid Mech. 9: 193-217. 


. Phillips ,O. M. (1978) The Dynamics of the Upper Ocean, 
Cambridge University Press. 


. Phillips, O.M. (1981). J. Fluid Mech. 107: 465-485. 


. Pokazeyeu, K.V. & Rosenberg, A.D. (1983). Oceanolo- 
giya 123: 429-435. 


. Rockliff, N. (1978). Math. Proc. Cambridge Phil. Soc. 83: 
463-479. 


. Shemdin, O.H. (1990). J. Geophys. Res. 95: c9. 16319- 
16338. 


. Shyu, J.H. & Phillips, O.M. (1992). J. Fluid Mech. 217: 
115-141. 


Naval Research Reviews 


. Su, C.L. (1973). J. Waterways Harbor Coastal Eng. Div. 


ASCE 375-391. 


. Thews, J.G. & Landweber, L. (1935). U.S. Expr. Model 


Basin Report 408, Navy Yard, Washington, D.C. 


. Trulsen, K. & Dysthe, K.B. (1990). Water Wave Kinemat- 


ics A.Torum, O1.,T. Gudmestad (eds.) Kluwer Academic. 


. Trulsen, K. & Dysthe, K.B. (1992). in Breaking Waves, 


IUTAM Symposium in Sydney/Australia 1991, M.L. 
Banner, R.H.J. Grimshaw (Eds.), Springer-Verlag Berlin 
Heidelberg. 


. Trulsen, K. & Mei, C.C. (1993). J. Fluid Mech. 251: 


239-271. 


. Trulsen, K. & Mei, C.C. (1995). J. Fluid Mech. 290: 


345-376. 


. Tulin, M. & Miloh, T. (1990). Proc. 18th Symp. Naval 


Hydrod., Ann Arbor. Ursell, F. (1953). Proc Cambridge 
Philos. Soc. 49: 685-694. 


. Wilton, J.R. (1915). Phil. Mag. 29: 688-700. 
. Wa, TY. (1987). J. Fluid Mech. 184: 75-99. 
. Wu, J.-K. & Liu, P-L.-F. (1990). J. Fluid Mech. 217: 


595-613. 


. Wu, D.M. & Wu, T.Y. (1982). Proc. 14th Symp. Naval 


Hydrod., Ann Arbor, 103-129. 


. Yuen, H.C. & Lake, B.M. (1980). Ann. Rev. Fluid Mech. 


12: 303-334. 


. Zabusky, N.J. & Kruskal, M.D. (1965). Phy. Rev. Lett. 15: 


240. 


. Zhang, J. & Melville, W.K. (1990). J. Fluid Mech. 214: 


321-346. 


. Zhou, C. & Liu, P.L.-F. (1987). J. Fluid Mech. 175: 


143-155. 





What Can We Learn from 
Wave Tank Experiments? 


Norden E. Huang and Steven R. Long 
NASA Goddard Space Flight Center 


Laboratory for Hydrospheric Processes | Ocean and Ice Branch 


Greenbelt, MD 


Abstract 


Wave tank experiments have provided the indispensable 
data for the past advances in the studies of nonlinear dynamics 
of the ocean waves. Experimental results on nonlinear wave- 
wave interactions and wave instabilities, wave breakings, wave 
Statistics, wave induced longshore currents, wave blockings 
and trapping when encounter currents, and wave observation 
under nonstationary conditions all showed the importance of 
the laboratory observations. In most of the cases, field obser- 
vation just cannot provide detailed enough data for us to 
understand the physics of the phenomena. Thus, laboratory 
experiments are not only indispensable, but also the most 
cost-effective way to study the nonlinear ocean wave phenom- 
ena. 


Introduction 


Wave motion was among the first natural phenomena 
studied mathematically by scientists. The giants in mathemat- 
ics and theoretical mechanics such as Lagrange, Airy, Stokes, 
Laplace, Russel, and Rayleigh have all left their marks in this 
field. Past studies, however, have been centered on the infini- 
tesimal undulation of the water surface with harmonic distor- 
tions. New and decper understanding of wave dynamics only 
started with the discoveries of weak nonlinear wave-wave 


interactions that describe the longer term wave evolution. With 
few exceptions, most of the theoretical advancements made 
were inspired by direction observations with the laboratory 
observations constituting most of the success. These observa- 
tions provided not only confirmations of the existing theories, 
but also new insights and inspirations for further develop- 
ments. The field observations, though indispensable for final 
confirmation of the theory for applications, will probably 
never play as crucial a role as the laboratory in advancing the 
fundamental understanding of wave dynamics. 

Being motions confined to the air-water interface with a 
huge density difference between the air and the water, even the 
highly nonlinear wave motions are amenable to the inviscid 
approximation. Consequently, an irrotational motion can offer 
an accurate depiction of the phenomenon. Because of these 
properties, the laboratory experiments can provide good simu- 
lations for us to study most of the controlling mechanisms of 
the wave dynamics. Nevertheless, as close as the laboratory 
results can be to those of the natural phenomena, a laboratory 
experiment still cannot fully simulate the complicated balance 
of different natural forces. Therefore, in designing an experi- 
ment, one should know what mechanism to simulate, what 
forces need to be controlled and monitored, before mounting 
the tests. Through careful planning, detailed measurements 
under controlled conditions can yield enough data for compari- 
sons with the theory, and for reaching statistically significant 
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Figure 1. 


Resonant wave interaction measured by McGoldrick et al. 
(1966). Wave with frequencies f, and f, were generated 
mechanically. According to Phillips (1960) theory, resonance 
will occur when f /f,=1.775, at the frequency 2f,-f,, The two 
spectra measured in a wave tank showed the resonant 
interaction at work. 








a (mmm) 


























a (mm) 



































f (eycles/sec) 





conclusions. Additionally, laboratory experiments are ex- 
tremely attractive in the present limited budgetary period, for 
an experiment can be completed with a minimum amount of 
expense comparing to the full field tests. In the following 
sections, we will give a few examples of the wave properties 
we have learned from the laboratory experiments, and the 
insight gained from such observations. 


Laboratory Experiments 


Nonlinear Wave-Wave Interactions and 
Instability 


The starting point of modern water wave theory can be 
placed rightfully at the discovery of the weakly resonant 
nonlinear wave-wave interactions by Phillips (1960). Later, 
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Lighthill ( 1965, 1967) pointed out that a deep water wave train 
can be unstable. Phillips (1967) identified the instability as a 
consequence of the nonlinear resonant interactions. The non- 
linear interaction and the wave instability thus became the 
central problem for the deep water wave studies. Laboratory 
experiments by McGoldrick, et al. (1966) and Longuet-Hig- 
gins and Smith (1966) provided quantitative evidence of these 
interactions as shown in Figure 1. The theoretical and experi- 
mental observations by Benjamin and Feir (1967) and Ben- 
jamin (1967) have provided conclusive proof of the instability, 
as well as demonstrate how theory and experiment can coop- 
erate to advance our knowledge into new territory. Indeed most 
of the salient properties of the wave evolution such as the 
instability and the frequency down shift (Hasselmann, 1962) 
were all observed in laboratory settings (Lake et al. , 1977; 
Lake and Yuen, 1978; Yuen and Lake, 1982; and Huang et al., 
1986). These basic ideas are responsible for the nonlinear 
source function, the key part in the third generation wind-wave 
prediction model (Komen, et al, 1994). 

Following these pionecring studies, Su (1982 a, b), Su et 
al. (1982), and Su and Green (1984) have also observed higher 
order resonant interactions in the laboratory that involve five 
waves, resulting in the three-dimensional instability, as shown 
in Figure 2. These observational results have inspired the 
theoretical analysis of the higher order resonant interaction by 
McLean (1982 a, b). The theoretical and experimental obser- 
vations have defined the pattern of wave evolution to a high 
degree of accuracy. 

Even with these successes, the details of the frequency 
downshift is still unknown: Is it a continuous process as 
assumed in the classic treatment? Or is it a discrete process as 
observed by Lake and Yuen (1978) and Ramamonijiarisoa and 
Mollo-Christensen (1979)? The puzzle was further investi- 
gated by Huang (1995) and Huang et al. (1996), who proposed 





Figure 2. 


Three dimensional wave pattern generated at an outdoor wave 
channel by Su (1982). The initial waves are two dimensional. 














Figure 3. 


Wave elevation data measured in the laboratory with the wave 
maker set at 2.5 Hz, and an initial steepness at 0.2. The vertical 
axis indicates the distance from the wave maker. Notice the 
increasing modulation of the wave envelope as the waves 
evolve. 
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that the wave evolutionary processes were actually local, and 
abrupt. These conclusions were based on the analysis of data 
again from the laboratory by a pioneering new the Charac- 
teristic Scale Decomposition Method, which was then linked 
to their Hilbert Spectrum Analysis. This approach provides the 
local amplitudes and instantaneous frequencies that are finally 
used to produce a joint distribution that can be contoured as a 
time-frequency distribution. As shown by Huang et al. (1996), 
the Hilbert Spectrum produced a vastly better resolution in 
both time and frequency localization than the wavelet analysis. 

Let us examine the processes of frequency downshift in 
more details. Following Huang et al. (1996), who studied the 
phase changes starting from a uniform reference state. They 
showed that the wave frequency downshift was achieved 
through a local process: Two waves combine and become one, 
or three waves combine and become two. In general, n waves 
would combine to become (n - 1 ) waves. Other than the waves 
involved actively in these local processes, all the other waves 
are left unaffected. Each time, one and only one wave would 
disappear locally within the time span of only one wave period. 
The process is designed by Huang et al. (1996) as the ‘fusion.’ 
In the initial two-to-one stage, it is identical to the ‘lost crest, 
the ‘crest pairing’ and the ‘phase reversal’ idea reported before 
by others, but the subsequent observed variations were more 
general. They showed the process quantitatively through the 
phase function. 


The details of the wave fusion process are illustrated by 
the following laboratory data set. Figure 3 shows the original 
wave elevation data as used in Huang (1995) and Huang et al. 
(1996), with wave maker frequency of 2.5 Hz and the initial 
steepness set at 0.20. This gives an overall view of all the 
measurement locations, in which the vertical axis indicates the 
distance from the wave maker. The wave elevation data show 
that the waves at the shortest fetch (about 3 m) are quite 
uniform, but already have slight amplitude modulations. As the 
waves evolve, the amplitude modulations grow with the fetch. 
If one uses the Fourier spectral analysis to examine the fre- 
quency variations, he will only see the sideband growth, and, 
eventually, the lower lobe of the sideband overtakes the main 
peak to become the most energetic peak at station 8 (fora fetch 
at 8 m). 

The process of the frequency downshift was further ex- 
amined with the Hilbert Transform as in Huang et al. (1996). 
After applying the Hilbert transform to all the data, the un- 
wrapped phase is plotted in Figure 4. Because the overall slope 
of the phase function should be the mean frequency, the 
decrease in the overall slope of the phase function thus indi- 
cates a decrease of the mean frequency. At this scale, the trend 
seems to suggest a global, gradual and continuous change of 
the phase function or frequency. The slope of the first four 
stations, in fact, stays almost identical to that of the first station. 
The first time the phase function shows any deviation from this 
initially closely clustered group is station 5, but even this 
change in the phase seems only to be able to give a small and 
gradual change in frequency. 

Huang et al. (1996) examined the phase changes by taking 
the difference of the phase functions of all stations with respect 
to that of the first station as the reference state. The difference 
is shown in Figure 5. From this figure, we can sce that the 
changes of the phase functions consist of a series of steps with 
the sharp jumps confined in very short time spans. Thus, in 
seemingly stationary data, we find local transient phenomena. 
To have a quantitative assessment of the phase anomalies, they 
introduced the phase-amplitude diagram, a method new to 
water waves, which is a plot of the amplitude envelope as a 
function of the phase difference between the local values and 
that of the reference state. A typical result is shown in Figure 
6 for station 5. 

A section of the time series wave elevation data from 
station 5 is expanded and displayed in Figure 7 a and b in a 
two-way comparison. When the data from station 5 are super- 
imposed on its phase difference from station 1, the time 
domain associated with the jump in phase values could be 
easily confined within a single wave. Then the raw elevation 
data from station 5 are superimposed on those from station 1 
as shown in Figure 7 b. All the wave peaks obviously line up 
except at the location of the phase jump, where two waves fuse 
into one. Thus, in the time span of two waves seen at station 
1, now it contains only one wave locally at station 5. The period 
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Figure 4. 


Unwrap phase for all the stations as shown in Figure 3. At this 
scale, the change of phase values seems gradual and global 
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is doubled; the local frequency is, therefore, one half of the 
reference state. This local change in frequency also indicates 
group velocity change as reported by Melville (1983). These 
results have provided a quantitative measure to the qualitative 
observations by Lake and Yuen (1978), Ramamonjiarisoa and 
Mollo-Christensen (1979); they have also removed the uncer- 
tainty expressed by Melville (1983) on whether the period 
indeed doubled at the ‘crest pairing.’ The ‘crest pairing’ is 
clearly demonstrated as a combination of two waves, where 
the phase loss is precisely 2. 

Huang et al. (1996) have further concluded that the phase 
jumps shown here are cumulative, counting relative to the base 
of the reference state. On careful examination, one can see two 
types of phase jumps: (1.) The simple one involving two waves 
fusing into one; and the complex one involving n waves fusing 
into (n- 1 ) waves. The changes are local. Both characteristics 
are drastic departures from the traditional view of the wave 
motions: Other than the limited time duration where fusion is 
underway, the mean frequency remains unchanged from the 
initial state. The fusion region expands discretely and locally, 
and eventually covers the whole time axis. Each fusion event 
can be either a ‘crest pairing’ or a ‘trough pairing.’ The fre- 
quency downshift is an accumulated effect of fusions of two 
or more waves and a readjustment of local frequency. These 
processes are irreversible. This asymmetry is due to the pre- 
ferred direction of downward energy flow. These results also 
raise concerns about the past assumptions on slowly varying 
wave trains. These problems clearly need to be resolved theo- 
retically; we need a new paradigm for wave study. Wave tank 
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experiments have made a significant contribution to this so far, 
and the wave tank will continue to make significant contribu- 
tion in the future. 


Wave Breaking 


Wave breaking is the key in determining the wave energy 
dissipation, the turbulence generation in the upper ocean layer, 
the formation of the mixed layer, and the fluxes of energy and 
momentum across the air-sea interface. In fact, it can be linked 
to all the upper ocean layer dynamic phenomena as suggested 
by Huang (1986). Although great advances have been made in 
the last twenty years, progress in the study of breaking waves 
have been painfully slow. Our state-of-the-art understanding 
of breaking has recently been reviewed by Banner and Pere- 
grine (1993). It is a difficult problem to treat analytically and 
experimentally; the fundamental problems on breaking wave 
studies are still left unresolved: Of which the foremost is a 
precise definition of wave breaking. This seemingly obvious 
problem assumes a very different form in the wideband ran- 
dom sea, where breaking involves not only the obvious white- 
caps generating events, but also the countless not so obvious 
microscale breakers. Granted, it is widely accepted that a 
breaking event occurs when the velocity of a water particle 
near the surface exceeds the wave phase velocity. Still, there 
are so many velocity scales in a wideband wave field. The issue 
is further complicated by the fact that the surface velocity in 
the field is difficult to measure, while the only measurement 
of the surface velocity is again performed in the laboratory. 
Without a precise definition has hampered the measurements. 
As a result, optical detection of whitecaps, wave gauge detec- 





Figure 5. 


The phase function deviations from a reference state chosen 
as station 1. Notice the localized jumps of the phase functions 


at this scale. The amount of phase jumps are all 2x. 











Phase Difference (rad/(2*pi)) 


2000 3000 4000 


Time (sec) * 100 











Figure 6. 


The phases-amplitude diagram for station 5. Notice the dis- 
crete jumps in phase values. 
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tion of surface elevation jumps, radar detections of large Dop- 
pler frequency shift or large reflection spikes, and acoustic 
detection of noise generated by the bubbles generated by the 
breakers have all been used as surrogates for breaking with 
mixed results. Really quantitative measurements are actually 
taken place in the laboratory. 

Most of the laboratory breaking waves are generated 
mostly by the following ways: The first is in a steady flow with 
a submerged stationary hydrofoil as did by Banner and 
Melville (1976), Banner and Fooks (1985), and Banner (1990). 
A variation of this method is to tow a submerged hydrofoil in 
standing water as did by Duncan (1981, 1983). These studies 
established that the breaking wave can indeed increase radar 
reflectivity, generate undcrwater noise, enhance the momen- 
tum flux. The details of the breaking wave measurements also 
enabled Cointe (1987) to publish a theory of the dynamic of 
steady breaking waves, and established a numerical model. 
The second method is to generate breaking waves by focus 
waves in a preselected point as reported by Rapp and Melville 
(1990). Among the most significant findings is the great effi- 
ciency of breaking waves in transferring momentum and en- 
ergy from wave ficld to the water column. They have also 
measured the surface motion, momentum flux, energy loss due 
to breaking, breaking induced surface currents, turbulence 
generation, and mixing. These results constitute the most 
comprehensive body of our knowledge concerning the break- 
ing waves so far. 

On the theoretical side, Longuet-Higgins (1982, 1986), in 
a brilliant display of his analytic prowess, has found analytic 
solutions for overturning fluid motions. He showed that the tip 
of an overturning breaking waves can be model by a slender 


Figure 7. 


Detailed comparisons for a selected section of data at station 
5. 

(a) Raw data superimpoised on the corresponding phase 
function for the same time period. The time span covering the 
phase jump is confined within one wave period near the 49th 
second. 

(b) The same ra data superimposed on the data from station 
1 serving as a time scale here. The wave near the 49th second 
shows a clear 2-to-1 simple fusion event. 
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hyperbola. For a plunging breaker, the shape of the forward 
face can be modeled by a simple parametric cubic curve as in 
Figure 8. Such agreement between laboratory data and theo- 
retical results is indeed encouraging. Even with these remark- 
able successes, the difficult problem of understanding the true 
three dimensional breaking of the natural wind waves in field 
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The comparison between a laboratory plunging breaker front 
profile with the cubic parametric curve proposed by Longuet- 
Higgins (1986). 





density function of slope and elevation was also non-Gaussian, 
giving a further example of how theoretical development and 
laboratory verifications can work together to advance knowl- 
edge as in Figure 11. Encouraging as these earlier results were, 
they were still limited to narrow band cases, with building 
blocks of the theory resting on the Stokes wave. But the 
approach did give insight into the local statistical properties 
and the dominant mechanism, instantaneous modifications of 
the wave profile resulting in crest-trough asymmetry. Labora- 
tory data were also used to determine the key parameters in a 
wave spectral model as discussed by Huang et al (1990 a and 
b). 


Wave Induced Longshore Currents 


Longshore current is an important phenomenon in the 
littoral zone. Traditionally, the gencration mechanism has been 
identified as the breaking waves (see, for example, Longuet- 
Higgins, 1970 a and b). This mechanism is an important one. 
Neglected, however, is the role played by the non-breaking 
oblique long waves at the coastline, which has been shown by 
Zhang and Wu (1996) to be also important in generating the 
longshore current. Surface gravity waves, generated by wind 
either over the local region or at a great distance, can produce 
surf beats and drive the waves up on the beach. At an oblique 
incidence angle, the excess momentum can generate strong 
longshore currents even without breaking. To solve this prob- 
lem, one must consider the three-dimensional run-up of long 
waves on a beach with sloped bottom. The plane problems of 
run-up and reflection of waves incident normally on a sloping 
beach have been investigated by Carrier and Greenspan 





is still a formidable challenge. Our present knowledge con- 
ceming those events amounts only a theoretical model by 
Phillips (1985), and some detailed laboratory observation of 
the breaking wave kinematics by Bonmarin (1989). More 
quantitative data are needed for further advance our knowl- 
edge. 


Wave Statistics 


Wave statistics is important for engineering design, re- 
mote sensing, and for understanding the dynamics of the 
random wave field. Studies of the wave statistics are made 
difficult by the lack of data. Through theoretical analysis 
(Longuet-Higgins, 1963), it was found that the nonlinear effect 
can make the surface elevation distribution non-Gaussian. 
From laboratory data, Huang and Long (1980) demonstrated 
the agreement between theoretical result and actual measure- 
ments, as shown in Figure 9. This demonstrated that nonlinear 
analysis based on perturbation of potential motion of waves 
can work very well sometimes. This work was later extended 
by Huang et al. (1984) to demonstrate that the joint probability 
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Figure 9. 


Laboratory measurement of the probability density function for 
wind waves compared with the theory of Longuet-Higgins 
(1963) (smooth line) and the eight term Edgworth Expansion 
(dotted line). 




















Figure 10. 


Contour and prospective views of the joint probability density 
function for slope and elevation for a significant slope of 
0.0323. 
































(1958), Carrier (1966), and Tuck and Hwang (1972) based on 
Airy’s model. Numerical calculations have been tried by Kim 
et al. (1983) based on the Euler model. For the case of three- 
dimensional coastal dynamics, Carrier and Noiseux (1983) 
evaluated the oblique run-up and reflection of long waves on 
a plane beach based on linear wave theory. New results of 
oblique run-up, refraction and reflection over an arbitrary 
variable down slope were obtained by Zhang and Wu (1996) 
based on the same method as used by Teng and Wu (1992, 
1994) to the wave propagation in a channel of arbitrary shape. 
Most remarkably, the closed form solution of the wave induced 
mass transport near the coast without breaking has a magnitude 
comparable to that generated by the breaking ones especially 
at the water line. To sort out the details of the mechanism, a 
laboratory experiment was conducted at California Institute of 
Technology by them. A typical result is given in Figure 12 in 
which the dye streaks vividly demonstrated the strength of the 
circulation set up by the oblique incident waves. The magni- 
tude of the current agrees well with the theoretical predicted 
value. Unlike the breaking induced longshore current, this new 





Figure 11. 


Six consecutive top views of the longshore current generated by a 
train of obliquely incident cnoidal wave train with frequency at 1.23 Hz, 
at an angle of 29° to the coastal line. The beach is sloped at 7°. The 
six panels represent different time after the release of the dye: (1) 
t=0; (2) t=15 sec; (3) tf=28 sec; (4) t=45 sec; (5) t=66 sec; and (6) 
t=97 sec. The speed can be inferred by the movement of the dye 
streak over the bottom grid at 2X2 cm. 








model predicts the wave induced current to cover the whole 
nearshore zone. Recent ficld observation at Duck using tower- 
mounted video also indicated that the wave to current energy 
conversion process was not localized at the breaking point, but 
covered a considerable distance shoreward (Long and Sal- 
lenger, 1995). Thus, laboratory experiments again provided a 
demonstration of the mechanism for the important problem of 
long-shore current generation. This discovery is important for 
coastal circulation modeling, for it provides the boundary 
condition for the current field. In view of the importance 
assigned to the long ignored coastal zone dynamics, this phe- 
nomenon deserves careful further evaluation in the future. 


Wave Blocking and Trapping 


The blocking and trapping of wave packets by inhomo- 
geneous flow fields again have been studied in a laboratory 
setting by Lai et al (1989) and Long et al. (1993). Evidence of 
multiple reflections and thus trapping within the current gra- 
dient zone near the blocking point has also been reported by 
Pokazeyev and Rozenberg (1983). A strong current gradient 
zone near the blocking point was seen to trap wave energy, 
shifting the wave numbers into the capillary range as shown in 
Figure 12. The production of capillary waves by this mecha- 
nism is in the range of the microwave wavelength of many 
remote sensors, and can explain why SAR can image the 
bottom bathymetry, though the electromagnetic wave pulse of 
the remote sensor could not really penetrate the water surface 
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to any extent. When local current gradients are related to 
bottom features, then trapped waves can be found in patterns 
related to the bottom topography, and thus detected remotely. 


Nonstationary Wave Observations by 
Ensemble Averaging 


The laboratory is particularly well suited for the study of 
transient and nonstationary events. Under computer control, 
events such as the collision of two wave packets, the passage 
of a wave packet through a wind-wave field (see Chu et al., 
1992), and the initialization of surface roughness by sudden 
wind gusts (Spedding et al., 1992; Spedding, 1996) can be 
precisely repeated, which allows all the identical realizations 
necessary for good statistics to be studied by means of ensem- 
ble averaging techniques. As shown in Chu et al. (1992), this 
ensemble averaging of 100 realizations led to the discovery of 
some new phenomena, the overshoot of energy density as it 
recovers from the passage of a wave packet through a wind 
wave field, as shown in Figure 13. This overshoot was ob- 
served for all the conditions studied. In the surface roughness 
initialization study by Spedding (1996), the laboratory again 
provided the detailed quantitative data that would be almost 
impossible to collect in the field. Experiments like this can be 
conducted in a computer controlled laboratory as described by 
Long (1992). Such data can provide truly ensemble average of 
complicated nonstationary or transient phenomena. Thus, 
laboratory is posed to make another fundamental contribution 
in the air-sea interface dynamics. 


Conclusion 


All the above experiments have shown that wave tank 
experiments can indeed play a significant role in under- 
standing not only wave dynamics but also the associated 
dynamic phenomena. Such experiments can work with 
theoretical efforts to insure a real-world connection be- 
tween theoretical developments and measurable phenom- 
ena. It is the classical cause and effect relationship in action: 
theoretical efforts can examine the cause and predict the 
effect, while experimental methods can verify the predicted 
effect, and even suggest alternate causes. This can confirm 
theory, or reveal reasons for its revision or expansion along 
another path. To date, most significant wave related dy- 
namic phenomena were first found in the laboratory. There- 
fore, the laboratory studies are indispensable means for 
water wave investigations. 

Wave tank experiments are also extremely cost-effec- 
tive. Computer controlled facilities are available now for 
studying the most complicated nonstationary and transient 
wave phenomena. In the present tight budget situation, 
laboratory study should be seriously considered as an alter- 
native. 
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Figure 12. 


A typical case of non-breaking blocking of waves by opposing 
current. The waves propagates from left to right encounter a 
current of 30 cm/sec. The waves are stopped without 
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point as indicated by the fine wave pattern reflection on the 
back wall. 


Figure 13. 


The mean-square surface displacement of short wind-waves 
passing a fixed elevation probe. The is the ensemble average 
of 100 realizations, with the ensemble average of the longer 
wave packet removed. The wave packet passes this location 
between the 8 and 14 second mark 
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